function [pi] = profit_Sh(pi_g,pi_b,m,mu,p_g,p_b,lambda)

    p=mu*p_g+(1-mu)*p_b;
    np=1-p;
    gamma_g=mu*p_g/p;
    gamma_b=mu*(1-p_g)/np;
    pi=-1*(mu*p_g*(pi_g+m)+(1-mu)*p_b*(pi_b+m)-lambda*ShannonCost(mu,p,gamma_g,gamma_b));